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Abstract 



We establish a characterization of the class of the simplest nearest neighbor spin sys- 
tems possessing the mean coverage function (mcf) that obeys a second order differential 
equation, and derive explicit expressions for the mcf's of the above models. Based on 
these expressions, the problem of ergodicity of the models is studied and bounds for 
their spectral gaps are obtained. 



1 Introduction and Summary 

It is commonly acknowledged that even in the case of a simple infinitesimal interaction mech- 
anism, a description of the transient behavior of an interacting particle system (IPS) is an 
intractable mathematical problem in the theory of Markov processes. In view of this, there 
is a continuing interest, both in theory and applications, in seeking solvable (in some sense) 
models of IPS. One of the important functionals of IPS is undoubtedly, their mean coverage 
function. In the present paper we continue to study the behavior in time of the mean cover- 
age function of a class of IPS called the simplest nearest neighbor spin systems (SNNSS) on 
s- regular graphs. Namely, developing the approach of Granovsky and Rozov , we establish 
a characterization of the class of SNNSS that posses a mean coverage function satisfying a 
second order differential equation. This is the main result of the present paper, stated in the 
Theorem, Section 3. The theorem asserts that the above class consists of the following four 
different modifications of the basic voter model: noisy voter model, noisy voter model with 
threshold = 2 ( or 3) on 2- (resp.3) regular graphs, a special case of a general threshold = 2 
model in one dimension and a degenerate model with thresholds on s-regular graphs. 
It should be noted that the first of these models is the unique SNNSS that has a mean coverage 



function satisfying a first order differential equation. This was proven in [pL2| . In Section 3 
we derive explicit expressions for the mean coverage functions of the above four models, by 
solving the corresponding second order differential equations. The formulae obtained show 
that adding a constant noise to flip rates results in considerable change in transient behavior 
of the process. This matter is discussed in Section 5. 

The next two sections are based on the aforementioned formulae for the mean coverage func- 
tion. Section 4 is devoted to the mean density function. We prove here that, when started 
from the product Bernoulli measure, the mean density functions of the above processes do 
not depend on the size of the graph. This remarkable property is used for the study of ergod- 
icity in the next section. Section 4 contains also a historical sketch of research related to the 
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subject. 

It is clear that transient behavior of the mean coverage function, which is of interest in itself, 
also provides information on the long-time properties of the process considered. In view of this, 
the last Section 5 is devoted to ergodicity and bounding the spectral gap for the class of models 
defined in Theorem. We give here a positive answer on the open problem about ergodicity of 
threshold = 3 noisy voter model on 3-regular graphs, for some values of parameters. Based on 
the expressions derived in Section 3, we obtain the upper bounds for the spectral gap of the 
four SNNSS. These upper bounds are compared with the lower bounds given by the e — M > 
condition (for references see [H], p. 31). 

Finally, note that in the course of the proof of the Theorem we derived identities that hold 
for a coverage of sites of a regular graph by O's and l's. These identities might be helpful in 
the study of other problems related to time dynamics of SNNSS. 

Most of the notation and language of our paper have been adopted from the seminal mono- 
graph on IPS |14]], by Liggett. 



2 Background 

We consider throughout the paper a SNNSS on a s-regular graph G of finite size N, with the 
set of vertices (sites) V = {x}. Recall that a graph is called s-regular if each of its vertices has 
s neighbors. By SNNSS we mean a time homogeneous Markov process ip t , t > with state 
space = {0, 1} V = {77} and the infinitesimal time dynamics given by (|2.2| ) below. The 
elements r\ = {rj(x), x G V} of X^ are called configurations. We will say that a site x G V 
is occupied(resp., empty) in the configuration rj G X^, if rj(x) is 1 (resp., 0). The SNNSS are 
featured by the property that the flip rate c( x, rf) of a spin at a site x G V in a configuration 
rj G X^ depends only on the number k(x,rj) of occupied neighbors of x in the configuration 
77. Formally, 

c(x,T)) = A fe (l - T](x)) +/i k r)(x), k = k(x,7)), x G V, rj G X N , (2.1) 

where k = 0,1, ... ,s (resp., k — 0, 1, . . . , s ) are the rates of the infinitesimal 

transitions — > 1 (resp., 1 — > 0) at a given site in a given configuration. Finally, denoting by 
r) x the configuration obtained from rj by flipping the spin at the site x, the above assumptions 
conform to the following infinitesimal time dynamics of (p t , t > : 

Pr((f t+ At = i]x I ft = rj) = V)^t + o(At), At > 0, t > 0, xeV, rteX N , (2.2) 
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o(At) 
At 

So, to compare with a variety of the so called biased models, ( see e.g. Madras, Schinazi and 



where ~~ * ^> as ^ ~~ * ^- 



Schonmann |]16|) SNNSS is a spatially homogeneous process. 

It is known that in the above setting the process ip t , t > is fully defined by the 2s + 2 
parameters > 0, fik > 0, fc = 0, . . . , s. Namely, the generator Q of the process is given by 



fi/fa) = E c MX/fa*)-/fa)), /eC(AV), ^4 (2-3) 

where C(X N ) is the class of bounded functions / : X N — > R. 

Denote (ft , t > 0, 77 G the SNNSS starting from a configuration 77 and 



Mf\t) = Ef((p[ v) ), t > 0, / G C(^jv), 77 G ^v. As in Granovsky and Rozov p[, our starting 



7 

point will be the following assertion that is a straightforward consequence of the Hille- Yosida 
theorem. 

Proposition 1. The function Mf (t), £ > satisfies, for all 77 G Afjv, a linear differential 
equation of order /, (/ > 1) 

— ^ = 5^— ^ + 5, t>0, V eX N (2.4) 

with coefficients Ai, i = 0, . . . , I — 1 and B that do not depend on 77 G A'/v and t > 0, iff the 
generator Q of the Markov process considered obeys the condition 

Q l f = J2AiWf + B, 7]&X N , (2.5) 
where := £ = 0, 1, 

Our subsequent study of the characterization problem described in the previous section is 
based on the fact that (|2.4|) is equivalent to (|2.5|). 



3 Main result 

The coverage of the graph G by a configuration rj G Xn is the function \rj\ : Xn — > -R + , |t/| = 
Exev^W and M (r?) (t) := E\<pl v) \, t>0,r]eX N is called the mean coverage function of the 
process <fi\ t > 0. The function M^(t),t > is one of the most important functionals in 
applications. In Granovsky, Rolski, Woyczinski and Mann || and Belitsky, Granovsky || the 
function was studied in the context of adsorption - desorption process given by 
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Afc = A > 0, Hk > 0, k = 0, . . . , s. It was observed there that the function M^ v '(t),t > 0, has 
a saddle point, under certain conditions on parameters of the process. 

Our main objective will be to describe the class of SNNSS satisfying (|2.4j ) with / = \rj\ and 



I = 2. For I = 1 the problem was posed and solved in [E|. 
We introduce some more notation. Denote 

9i {rj) :=^(M), z = 0,l,..., ?7 e X N (3.6) 

to obtain from (|2.3|) 

</ifa) = E^X 1 " 2 ^)). (3.7) 
In view of our objective, we will need to unlock the structure of gi . 

Let D be a nonempty subset of V. We will say that y G V is a neighbor of D : y ~ D, ify^D 
and ?/ is a neighbor of at least one site in D, and we denote S\(D) the set of all neighbors of 
the subset D. In particular, by Si(x) we denote the neighborhood of x G V. We also define 
5i(x) = 5i(5j_i(x)), £ = 1, 2, . . . , 5q(x) = {x}, x G V. 
Next, for any x G V define the difference operator A x : C(Xjf) — > C(Xn) 

&*f(v) = f{Vx)-m, feC(x N ), xev (3.8) 



and write A® / = A x A y f, f e C(X N ), x, y G V. Then, by our definition (|3~6p and (O 



we have 



#2(7?) = E c(x,7])A x g 1 (7]), r] G ^jv. (3.9) 



Further, it follows from (13.71) and (12. II) that 



xe8 1 (y) 

(c(y, rj y ) + c(y, rj)) (1 - 2rj(y)), rj G X N , (3.10) 



for any y G V. Since A^ = A y ^ x , x,y G V, ( 3.10Q implies the important fact that 



Ag)^ 1 ( 77 ) = 0, 77 e**, (3.11) 
whenever x G" <5i (2/) U^d/) and x ^ y. This and ( |3.9| ) give 
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A y g 2 (v) = A y [c(y,r])A y g 1 (r})\ + 

E a 

y \c(x,T])A x gx(r])\, y G V, rj G X N . (3-12) 

For the proof of our main result, stated in the Theorem in the sequel, we need to impose the 
following two conditions on s-regular graphs G considered. 

(i.) First, we assume that G is triangular free graph, which means that if x, y, z G V : y, z ~ x, 
then y, z are not neighbors. The second condition is a technical one. 

(ii.) We assume the existence of a pair of vertices y, z G V s.t. z G 6 3 (y) and the two sets of 
vertices Ei t2 := Si(y) fl 5 2 (z) and E 2) i : = S 2 (y) H &l(z) are singletons. 

Observe that the conditions (i) and (ii) are satisfied e.g., when G is an s-regular tree or 
G = Z d , d > 1. 

For the purpose of establishing our characterization result we employ a technique that is 
presented below. We start with the notations adopted from || and |12| . Denote = 



n k(v)i k — 0, . . . , s, i — 0, 1 the number of occupied (i=l) (resp., empty (i=0)) sites having 
k occupied neighbors in a configuration r] G X n , and let = nf\r]), = n^^r/), k = 
0, . . . , s denote the corresponding sets of vertices x G V in a configuration r/ G A^r. Finally, 
we denote = V^(rf),i = 0, 1 the set of all empty (resp. occupied) sites in rj G Xn- 
Then gi defined by ( pTT|) can be expressed as 

s 

9i(v) = E( A fc4 0) - ton* 77 G <*jv- (3.13) 

fc=0 

The following identities that are valid for any s-regular graph will be crucial for our subsequent 
study: 

P = P( V ):= £ A x n« = £ A X 0) ; V^Xn (3-14) 



Qo = Qo(v)- = E A^ ) = -n(°)+ni 0) 1 , t? G ^ (3.15) 

Qi = QM-= E A x n^ = -S7if +n« 77 e** (3.16) 

xevW 

z ^A a nW) = -n?), £ A^nf* = -n? 3 , /? G (3.17) 



r- (°) r- t 1 ) 
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The proof of the identities can be obtained after some thought from the preceding definitions. 
We also write 



Ri = Ri(v) := E A * n o\ i = 0, 1, V e (3.18) 

Let 77 be the configuration obtained by flipping the spins at all sites x G V in a configura- 
tion G Afjv- Then we have rajj. (77) = n^~^\r]), k — 0, 1 . . . , s, i = 0, 1, 77 G Afjy and, 
consequently, 

Ri(v)=Qi-i(fj), r/GAV, « = 0,1. (3.19) 

Now we are in a position to state the following 
Lemma. The identity (in 77) 

n ( P + nf\ - - nf = (nf - nf^Fi + \v\F2 + F 3 , rj G (3.20) 

where Fj, i — 1, 2, 3 are coefficients that do not depend on 77, holds iff s = 2, 3. In both cases 
of s, F 2 = 2, F 3 = -TV, while 

' -3, if s = 2 , 

F x = { 3.21 
-2, if 5 = 3. 



Proof. We put in ( |3.20|) first 77 = and then 77 = to find F 2 and _F 3 . Now consider the 



case s > 3. Due to the fact that the graph considered is triangular free, we have for any 
x ~ y, x, y G V, 

n ( }\Kv) = nflAKy) = 4\Kv) = nf\Kv) = 0, 4>^) = N - 2s, n?(V) X;y ) = 2. 

(3.22) 

Substituting this in ( ^20|) , gives -2 - (JV- 2s) = 4- AT, which says that ( ^20|) does not hold 



for s > 3. 

Further, if s = 3 and 77 = $ x , x G V, then 

4 1) (0x)=4° ) (0x) = 4 1) (0 :c )=4 O) (0x) = O, 7i^ o) (0 x ) = AT-4, 4 1) (0,) = 1. (3.23) 
The latter implies F\ = —2. 
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In the case s = 2 we have n^°_ 1 (0 Xi2/ ) = 2, and the same argument as before gives F\ = —3. 
Finally, it is left to show that the identity Q3.20| ) indeed holds for s = 2, 3. We use the 
relationship 



2^1 - N = \rj\ -(N-\ V \) = - 4 0) ; 

fc=0 



to obtain for s = 3 



n 



+n W_ n W_ n W + 2( n W- n W)-2\r)\+N = ^2{knl l) -{s-k)nl l> ) = 0, V G ^, (3.24) 

fc=0 

where the last equation follows from the identity J2t=i k(nf^ + ) = s\r]\, r] G X n that is 
valid for all s-regular graphs. The same argument proves the assertion for s = 2. £ 

Finally, we will distinguish the following modifications of the Basic Voter model : 

Noisy Voter Model. Afc — Afc_i = fik-i ~ — d, k = 1, . . . , s. The model was introduced in 



12l and intensively studied in [ llfl . Here the noise is given by the two parameters hi = Aq, h 2 



/Jo — sd added to the basic voter model (see p^|,[^q|): A^ = kd, \ik = (s — k)d, k = 0, 1 



Note that in |TJ| Ex.2. 5, p. 136, it is considered a general (i.e. not necessarily the nearest 

neighbor) version of voter model with noise. 

Noisy Voter Model with Threshold = q (1 < q < s). 



A fc = /ifc+s-g+i = h> 0, k = 0, . . . ,q - 1, \ k = /j fc _ 9 = h + a > 0, k = q,...,s. (3.25) 

This is the simplest case of a nonlinear voter model. In the case h = (the absence of noise), 
the model was suggested by Cox and Durrett in J7J. ( For updated references see |TB|). In J7| 
it was also considered the threshold voter model with noise added to the death rates only. If 
q = s, then by scaling all the rates by the factor (2h + a)" 1 the model becomes the nearest 
neighbor Majority Vote Process (pjl, Ex.4.3(e),p.33 and Ex. 2.12,p. 140). 
Generalized Threshold Model with threshold — q (1 < q < s). The model is obtained 
from the previous one by adding a constant either to s — q + 1 birth rates A&, k = q, . . . , s, 
or to s — q + 1 death rates k — 0, . . . , q — s. Explicitly, 

Afc = fik+s-q+i = h > 0, k = Q,...,q-l, Afc = h + a > 0, k = q, . . . , s, 

iik = h + b > 0, k = 0, . . . , s - q. (3.26) 
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Theorem 

The mean coverage function M^"(t), t > of a SNNSS ip t , t > satisfies, for all 77 G X N , a 
second order linear differential equation 



A 1 TT ^ + A M {ri \t) + B, t>0, 7]eX N 



(3.27) 



dt 2 " x eft 

with coefficients A , Ai, B that do not depend on r\ G Afjy and t > 0, iff <^ t , t > is one of the 
following four models (Ci) — (C4) : 
(Ci) A noisy voter model. 

(C2) A generalized threshold model with threshold = s and h = ab = or h = 0, a = b. 
(C3) A threshold noisy voter model with threshold = s, when s = 2, 3 

(C4) A generalized threshold model with threshold s = 2 and h,a,b : h(a + b) = ab, h > 
0, h + a > 0, + 6 > 0. 

Proof. By virtue of Proposition 1, (|3.27 ) is equivalent to 



g 2 (< n )=A 1 g 1 {r ] )+A \r}\ + B, rj G X 



N ■, 



(3.28) 



The main difficulty is to prove that ( |3.28j ) implies one of the four conditions (Ci) — (C4) on 

the rates of <pt, t > 

If (ggp holds, then, by (gTTD 



A^g 2 ( V ) = A 1 A y %(r ] )=0, y £ V, z e 5 3 (y), V G X 
From the other hand, we get from (ft.l2|) and ( p.ll| ) 



& y %(v)= E A y 2 l[c{x,r ] )A x g 1 { V )}, y G V, z G 6 3 {y), V e X N . 



In view of ( |3.11| ) this gives 



(3.29) 



(3.30) 



& y %(v)= E ^y[c(x,v)^9M] + E ^[c(x,r,)A^ gi (r])], ze6 3 (y), rj G X N , 

(3.31) 

where we denoted : = Si(y) fl (^(z) and #2,1 := S 2 (y) H £1(2). We also derive from (|3.10|) 



Agki(»7)= E (l-27?(«))Ag) c ( w>77 ), xei 2 (z), r/G^ (3.32) 
ue5i(x)nSi(z) 



and 
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^ y 9i(v)= E (l-2r ] (u))A^ y c(u, V ), xe6 2 (y), rj G X N . (3.33) 

ueS 1 (x)n8 1 (i J ) 

We substitute now these expressions in ( |3.31| ) to obtain 



^MV) = E *v[<*,ri) E (1 - 2r ? ( M ))Ai 2 ) C ( M ,r ? )] + 
E A z [c(x,r/) £ (l-2r / ( M ))A( 2 ) C ( M ,r ? )], J/GV, 2G%), qe4 (3.34) 

xGi?2,i «£5i(a:)n5i(i/) 

Our immediate aim is to find conditions on the parameters of a SNNSS, imposed by the 
requirement 

Agl&fa) = o, )/ey, zeisd/), ije4 (3.35) 

Let, in accordance with the assumption (ii) on G, the vertices y, z in ( |3.34| ) be such that 



E 1>2 = { Ul }, E 2il = {u 2 }, (3.36) 
where Ui,u 2 G V. Then, in view of the above definition of the vertices Ui,u 2 , ( |3.34| ) becomes 



A^g 2 ( V ) = (A y c( Ul , V ))(l-2 V (u 2 ))A^ z c^2,v) 

+ (A z c{u 2 ,r))){l-2r ] {u 1 ))AW y c{u 1 ,r ] ), 77 G X N , ze5 3 (y). (3.37) 

The last expression will be our main tool in the subsequent study. 
We see from (|3~37p that 

Ag<? 2 (%) = -AgHpifo), A^g 2 ( Vz ) = -A^g 2 ( V ), z G 5 3 (y), 77 G X N . 



)k+l 



In view of this we set in (|3.37|) , 77(7/) = r](z) = . We also agree to write A(«)& = 
and A (2) (») fc = (»)ah-2 - 2(m) k+1 + (•) Jfe) where (•) is either A or fi. 

Now we will be attempting to find the explicit form of ( |3.37| ) in the following three cases of 
77 G Xfq that exhaust all the possibilities. For brevity, we denote fcj = k(ui,rj), 77 G X^. It is 
important to note that 8\(u\) 0^1(^2) is the empty set, since G is triangular free. 
Case 1. 77 G X N : 77(111) = 77(1*2) = 0. 



A(% 2 (r/) = (AA fcl )(A( 2 )A fc2 ) + (AA fc2 )(A( 2 )A fcl ), zeS 3 (y), 77 G X N , 



(3.38) 
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where < ki, k 2 < s — 2. 

Case 2. 77 E X N : 77(1*1) = 77(^2) = 1- 

Agi&fa) = (A^A^-i) + (A/i fe2 )(A( 2 V fcl -i), z e 5 3 (y), 77 G (3.39) 

where 1 < k\ , k 2 < s — 1 . 

Case 3. rj E X N : 77(111) = 0, 77(1*2) = 1. 

Agifefo) = (AA fcl )(A( 2 )/i fc2 ) + (A Mfc9 )(A«A fe _ 1 ), (3.40) 
where 1 < kx < s - 1, < k 2 < s - 2. 

We now know from the three cases considered, that the condition (|3.35|) implies 



(AA fcl )(A( 2 )A fe2 ) + (AA fc2 )(A< 2 )A fel ) =0, < h, k 2 < s - 2 (3.41) 

(A/i fcl )(A( 2 V fc2 -i) + (A^ fc2 )(A( 2 V fcl -i) =0, l<h,k 2 <s-l (3.42) 
(AA fcl )(A( 2 ) / u fe2 ) + (A /Ufc2 )(A( 2 )A fcl _ 1 ) =0, 1 < h < s - 1, 

0<k 2 <s-2. (3.43) 



Setting in ( ggg ) k x = k 2 = k gives 

(AA fe )(A (2) Afc) =0, 0<fc<s-2 (3.44) 
(A/i fc ) (A( 2 Vfc-i) =0, 1 < A: < s - 1. (3.45) 



Using this fact, we multiply the equations (p.41|) , (|3~4"2]) by AA^ and by A\i kl correspondingly, 
to obtain 



AA fel )(A( 2 'A fe2 ) =0, 0<k 1 ,k 2 <s-2 
A/i fcl )(A (2 Wi) =0, \<k 1} k 2 <s-l. (3.46) 



Since, by our definition, A^ 2 ^(«) fe2 = A(«) fc2+1 — A(») fc2 , the equations Q3.46| ) imply for s > 2 

A (2) A fc = 0, 0<A;<s-3. 

A (2 Vfc = 0, 1 < k < s - 2. (3.47) 
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Thus, we have for s > 2 



AA fe : 
A/x fc : 



d\, k — 0, . . . , s — 2 



d/j.1 k — 1, . . . , s 1 



Finally, in view of ( $A7\) and QOSD , we obtain from (g4T|)-(g4|) for all s > 2 



(3.48) 



d x A^X 



s-2 



^A (2 Vo 








(AA s _ 1 )(a( 2 Vo+(A/x )(a (2) A s _ 2 ) = 0. 



(3.49) 



Summarizing the preceding argument we conclude that ( |3.48|) together with Q3.49Q are nec- 
essary and sufficient for ( |3.35|) . Our next step will be devoted to show that the conditions 
Q3.48j ), (|3.49|) on the parameters A&, /i*,, k — 0, . . . , s imply one of the conditions (Ci) — (C4). 
Assume first that in fl3.49| ) d\ 7^ 0. Then we should have A( 2 )A S _ 2 = AA s _i — AA S _ 2 = 0, and, 
consequently, in view of ( |3.48| ), AA s _i = d\ 7^ 0. Hence, in view of the last equation in ( p.49|) , 
we obtain 

A( 2 Vo = A( 2 )A s _ 2 = 0. (3.50) 

By the same argument, (|3.50|) should also hold under the assumption d^ 7^ 0. ( |3.50|) together 
with (|3.48| ) is equivalent to saying that the flip rates A&, k — 0, . . . , s and /i^, k — 0, . . . , s 
form arithmetical progressions: 



AA fc = d x , k = 0,...,s-l 

A/i fc = d^ k = 0,...,s-l. (3.51) 

So, assuming d\ 7^ 0, one has 

A y c(x,T]) = (1 -2r)(y))(d x (l -t](x)) + d ti r)(x)), xeS^y), x,yeV, rj e X N . (3.52) 
Now (|3~52|) and ( gTOl) imply 

A y g 1 ( v ) = (l-2r)(y))(sd x -2k(y, V )(d tl + d x )-\ -iio), V^V, 77 e X N , (3.53) 



A®c{x, 77) =0, y 2 e6 2 {y), xeV, rj e X J 



N; 



(3.54) 
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and 



^SMv) = 0, y 2 e 5 2 (y), 77 e X, 



AT- 



Hence, by virtue of (|3.28| ) it follows from (|3.55|) that 

A ?ly92(v) = 0, y 2 € 5 2 (y), 77 6 Afjy. 
By ( |3.51 ) the latter is equivalent to 



(3.55) 



(3.56) 



^SMvi^yiQM) = -4(d A + <y (1 - 2t 7 ( 2/ ))(1 - 2 V (y 2 )) ((1 - 77^))^ - TjfeiK) = 0, 

1/2 e 5 2 (y), yi G <Ji(y), 77 e <*V- (3-57) 
This implies d\ + d ll = 0, which by (|3.51| ), corresponds to the noisy voter model (Ci). Since 



the same conclusion is valid under the assumption 7^ 0, it is left to assume that d\ = = 0. 
In this case it follows from fl3.4S| ) that the parameters of the SNNSS are of the form 



A fc = A, k = 0, . . . , s - 1, /ifc = /i, k = l,...,s, X s = A + a, /j = // + &, (3.58) 
where a, 6 G i? are such that A + a > 0, + 6 > 0. 

Note that in the case considered all three conditions ( 3.49 ) are satisfied, because Q3.58 ) implies 



A^ 2 ^A S _ 2 = a, A^/Jq = —6. In view of the relationships Zlfc=o n i°' ) = N — \rj\ and 



EjUo^ = M> (CT) yields for the model ( fUgp 



#1 



(t;) = AiV - (A + /i)|r7| + anf - fcn^, 77 G 



A'- 



By ( |2.3|) we also have 



(3.59) 



/i 51 ^n^, 77 G AW, fc = 0,...,s, i = 0,l. 



sen, 



(1) 



(3.60) 



Next, we apply ( |3.60|) for nf^ and 7i[ ) 1 ' ) to obtain, with the help of the identities (|3.14|) - ( |3.17| ), 

n(nf ] ) = XQo + /1Q1 - anf +bP, r\ G X N , (3.61) 
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and 



Q(4 



0-h 



XR + uRi — bn$ + aP, r] G X N . 



Now we derive from ( p.59| ) the expression for g 2 that we will be working with: 



(3.62) 



92{v) = ~ (A + A i )fl , i( ? ?) + XaQ — fibRi + fxaQi — XbR — a 2 n^ + b 2 n§\ rj G X N . (3.63) 
Q3.63j ) and (|3.59|) show that for the model (|3.58|) the relationship ( |3.28|) holds iff 



T(rj) := XaQo - ^bRi + /iaQi - XbR - a 2 n^ + b 2 n { Q> - 



(an<°> - bn$>)A 2 - B 2 \rj\ - C 2 = 0, r/ G AW, 



(3.64) 



where A 2 ,B 2 , C 2 are coefficients that do not depend on r\ G Ajv. We put in (|3.64|) first 77 = 
and then 77 = to obtain C2 = —XbN and _B 2 = Ao + //a. 

We will treat separately the case s > 3 and the case s = 2. Since a = b = leads to a 
particular case of noisy voter model, we suppose in the sequel that a 2 + b 2 7^ 0. Let some fixed 
y, z G V obey the condition (ii), and ui, u 2 G V are defined as in ( [3.361) . 
The case s > 3. Consider the following two configurations: rji, defined by rji(y) = rji(ui) = 
0, 771(f) = 1, for all v 7^ Ui, y, and r) 2 = (771 )« 2 . It is eas Y to figure out the following relationships 



Qofa) = 2, i = l,2, Q 1 {r ]1 ) = N-2s, Q 1 (r f2 ) 
Ro(Vi) = RM = 0, n^\r li )=4\r ll ) 



N + 1 - 3s, 



0, i = 1,2. 



Substituting ( |3.65|) in ( 3.64 ) gives 



(3.65) 



and 



2Xa + fxa(N - 2s) - (N - 2)(A6 + //a) + XbN = 



2Aa + fia(N + 1 - 3s) - (N - 3)(A6 + //a) + AWV = 0, 



which implies 



(3.66) 



(3.67) 



fia(s — 1) = (a + o)A, a/x(3s — 4) = A(36 + 2a). 
By the same argument, applied to the configurations 771, 772 we also get 



(3.68) 
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\b(s - 1) = (a + b)n, A6(3s - 4) = //(3a + 26). 



(3.69) 



We will find all solutions of fl3.68| ) and ( |3.69D . First we see that A/i ^ implies a = b and 
consequently, A = /i > 0, s = 3. This gives the threshold voter model (C 3 ). If A/i = 0, then 
we should have Xa = Xb = fia = fib = 0. By ( |3.64j ) this implies a = b or ab = 0. In the first 
case, we have A = fi = 0, which is again (C3), while in the second case, A = fi = ab = 0, 
which is (C 2 ). 

The case s = 2. Taking r/i as above, gives 



Qofoi)=2, QifaHiV-'l, flifaH 2, i2o(»7i)=T4 O, (?7i)=n^(»7i) = 0. (3.70) 
Consequently, ( |3.64j ) implies (a + 6) (A — /i) = 0. If a + b = 0, then ( |3.64j ) becomes 



Aa(Q + #0) + MQi + #1) " a 2 (4° } - 4 1} 
(n 2 0) + nf^a - - A)M - AaiV = 0, i] E X J 



N- 



(3.71) 



Since in the case s = 2 



Q + i? = -2nf - 2< ; + N 



,(1) 



A/ 



and 



Qi + #1 = -24 0) - 2< ; + |r/|, 7] E X J 



(i) 



A'- 



we see that ( |3.71| ) implies a = 0, and, consequently, 6 = 0. 

Let now A = /i. Then, we employ (|3.18|) , (|3.15|) and (|3.16|) to rewrite (p.64|) as 



(3.72) 
(3.73) 



T(rj) = X[aL(rj) - bL(f})] - a 2 nf + b 2 n£ ] 
(an 2 0) - brff)A<2 = 0, 77 e Afr, 

where we denoted L(r/) = nj°^ + n 2 — | ^7 1 and A 2 = A 2 + 3A. So 



(3.74) 



T(T7) + T(fj) = (a - b) [\(L(r,) + L(fj)) - (a + b)(n { 2 0) + n^) 

(nf+nJftiJ =0, V EX N . 



(3.75) 
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First observe that if a = b and A = \i then we have the model (C3). Next, substituting in 
L{rj) + L(fj) = nf ] + n { 2 ] + + - N = -n { 2 ] - 



(Em 



we have 

(nf + n^ ] ){A 2 + a + b + A) = 0, 77 G X N . 
Consequently, A 2 = —a — b — A, which in view of ( |3.74|) yields 



(3.76) 



Aa^+ri^+nf-^D-Aft^+n^+n^-iV+^O+afe^-n^) = 0, 77 G (3.77) 



(o) „(i)^ 



A specific feature of the case s = 2 is that the following identity holds 



2( 



— 71q ) = rii — Hi , V 



So, we obtain from (|3.77| ) 



n ( 2 0) - n { o ] )(-Xa - Xb + ab) = 0, G ^ 



AT, 



(3.78) 



(3.79) 



which gives the model (C4). 

This completes the proof of the necessity of the conditions (Ci) — (C4). 

The proof that each of the conditions (Ci) — (C 4 ) is sufficient for Q3.27| ) is now simple. In the 



case (Ci) it was shown in [12] that 



gi{T}) = \ N - {\ + ti s )\v\, V^Xj 



V: 



(3.80) 



which implies g 2 (r)) = -(A + /i s )gi(r]), r) G X N . 

In the case (C 3 ) we have s = 2, 3, A& = fik+i — h, k = 0, . . . , s — 1, \ s = fi := h + a, 
where a E R : h + a > 0. So, ( |3.63|) becomes 



g 2 (r)) = -2hgM + ah(Q -R 1 + Q 1 - R ) - a 2 (n^ - rtf), 77 G ^ 



N ■ 



(3.81) 



By the Lemma, ( |3.15| ), ( |3.16| ) and ( |3 . 1 8| ) we further obtain for s = 2, 3 



02(77) = -2hgM + ha((F 1 ~s-l)(n^~4 ) ) + 2\r ] \-N)-a 2 (n^-n^ ) ), 77 G X N , (3.82) 
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where Fx is given by ( |3.21| ). Hence, in view of ( |3.59j ), we have in the both cases of s 

g 2 ( v ) = -(a + 8h)g 1 -6h 2 (2\r ] \-N), r) G X N . (3.83) 



Let now (C 4 ) hold. With the help of (|3.78|) it is easy to verify that (|3.74|) indeed holds with 
A2 = —a — b—h. Finally, in the case of the model (C2) we have either g2{v) — —bgi(r]), 77 e Xn 
or g 2 (rj) = -ag^rj), r] e X N . X 

Corollary 1. The mean coverage functions M^(t), r] e X N , t > of the models (Q), 

i — 1, 2, 3, 4 are given by the expressions (Dj), i = 1, 2, 3, 4 correspondingly: 

(Di): 



M^(t) = (H- _^_)exp(-(A + /i a )t) + T-5— , f>0, rj e X N (3.84) 

:d 2 ) : 

M^(t) = -nf)(r/) exp (-at) + \ V \ + n<® (77), t>0, rj & X N (3.85) 

M^(t) = n { o\ri) exp (-bt) + - njj 1 ^), t>0, T] e X N . (3.86) 



or 



P.): 



where 



M^(t) = de" 01 ' + C 2 e _a2 ' + — , t>0, r)eX N , (3.87) 



and 



:d 4 ): 



8h + a - (-l)V (8/1 + a) 2 - 48/i 2 
«i = ^ , i = 1, 2. (3.88) 

Ci = , C 2 = \T]\ - — - Ci, 7] e X N . (3.89) 

012 — ol\ l 



M™(t) = C ie -^ + C 2 e- a2t + N j k + a . t>0, V eX N , (3.90) 

on + a + 



where 



6A + a + 6 - (-l)\/4/i 2 + (a + 6) 2 + 8/i(a + b) 
Oi = ^ — , i = 1,2, ab = h(a + b) (3.91) 



and 



a 2 — "1 8/1 + a + 
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Proof. 

The assertions follow from the relationships (Ei), % = 1.2,3,4 below that hold for the models 

(Q), i— 1,2,3,4 correspondingly. 

(Ei): 

gM = \ N - (\ + {i s )\ril t)EX n (3.93) 



;e 2 ): 



or 



(E 3 ): 



g 2 {rj) = -agi(r)), a > 0, r) e X N (3.94) 

02(77) = -bg x (77), 6>0, V eX N (3.95) 
92 ( r] ) = -(a + 8h)g 1 -6h 2 (2\r ] \-N), V E X N (3.96) 



g 2 ( v ) = -(Qh + a + b) gi (r]) -h(8h + a + b)\ri\ + hN(4:h + a), 77 e X N . (3.97) 

Namely, the expressions (D!)-(D 4 ) are obtained by solving the second order differential equa- 
tions corresponding to (Ei)- (E4), under the initial conditions 

MW(0) = M, dM 2 (0) =9i(v), VeX N . (3.98) 

4 The mean density function 

Let v be a probability measure on the state space X N . Denote by ^\ t > the SNNSS 
starting from v ( this means that the distribution of is v), and denote by 
M^(t) = E u M^(t), t > the corresponding mean coverage function. The function w^\t) = 
N- l M^\t), t > is called the mean density coverage function corresponding to the initial 
distribution v. 

Historical remark. The mean density function was studied in a number of papers. In 
addition to the previously mentioned literature that is immediately related to the context 
of the present paper, we outline now some adjacent topics of research. Special attention 
was devoted to the contact process. Gray |13| investigated the behavior of the population 
profile function pt{x) := P((pi n \x)) — 1, x G V, t > 0, when G = Z and 77 is the empty 
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configuration flipped at the vertex x = 0. Belitsky [[[] treated a special case of the previously 
mentioned adsorption-desorption process, when /i > 0, = 0, k = 1,2 and G = Z. It was 
proven in that the function t > 0, where z^o is the measure concentrated on the 

empty configuration, possesses a saddle point. This extends the result of |J. Continuing the 
discussion in ||, Belitsky [1]] relates the above phenomenon to the violation of the classical 
Langmuir law, known in physical chemistry. Note, that from Jj]], as well as 0, one can see 
how complicated is the structure of the iterations i > 1 of the generator of the process 
considered. This explains the difficulties in the study of the transient behavior of functionals of 
contact process even in the case G = Z. The problem becomes much simpler in the framework 
of the mean-field theory, that corresponds to the case when G is a complete graph. In this 
case, a SNNSS conforms to the birth-death process (see Granovsky and Zeifman [|l(|). The 
limiting behavior, of the density process TV 1 1 | , as N — >• oo, when <p t , t > is the basic 
contact process, was extensively studied in the literature. For the most recent review of the 
topic see Durrett ||. 

In conclusion, we mention two papers devoted to voter models. Cox || derived the limit of 
the density process for the basic voter model on the torus in Z d , under an appropriate time 



scaling. Mountford JHJ considered a class of one- dimensional multitype IPS, that are featured 



by the following property of its generator Q : 

sup \tlf n (<Pt)\ < const, (4.99) 



n,r) 



where f n (n) = ^2 xe z:\x\<nV( x )- He proved that under condition (|4.99|) the coverage process 



fnivt), t > is a martingale plus a term that is negligible as n — > oo. In this sense these 
models can be viewed as a generalization of the basic voter model. It should be noted that the 
condition (|4.99p fails for all SNNSS (Ci)-(C 4 ), except only the case of the basic voter model. 



An important particular case of v is the product Bernoulli measure v v , < p < 1, defined by 



Mv) = II P n{x \ l -p) 1 ^ =pW(i-p) N - M , v e x N . (4.ioo) 

For a given < p < 1, the measure v v corresponds to the initial distribution on X N , such that 
all spins are i.i.d. Bernoulli random variables. In view of this, M^\t) = NE V [ t P \x), \/x e 
V, n G Xn, and consequently, w^ p \t) = E(p[ up (x), Vx 6 V, n e X^. Hence, the mean density 
w^ p \t), t > defines the marginal distribution of the process (p[ at any site x G V at time 
t > 0. It turns out that the densities w% if), t > corresponding to SNNSS's (Ci) — (C 4 ) , 
have the following remarkable property. 
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Proposition 2. The mean density functions w^ p \t), t > 0, < p < 1 of models (Ci) — (C4) 

do not depend on N. 

Proof. 

It follows from (ggf ) that 



^(77) = + A gi„ 2 (r}), i = 3, . . . , r} e X N , (4.101) 

where, in view of (Ei)- (E 4 ), the coefficients Aj, v4 do not depend neither on 77 G nor iV. 
We deduce from ( |3.7| ) that 

E^gxirj) = iV£ ( ^ ) [c(x,r/)(l-2r/(x))], Vx e V. (4.102) 

It is clear that the expected value in the RHS of ( |4.102p does not depend on N, for any SNNSS. 
Next, in (|3.28|) the coefficient B = NBq, where, by (Ei)- (E4), the factor Bq does not depend 
on N. So, (|3.28|) implies E ( - Up ^g2(r]) = Nq2, where the factor does not depend on N. This 
together with (|4.101|) and (|4.102|) gives E^g^r/) = Nqi, i = 3, . . . , where again the factors 
q i: i = 3, . . . , do not depend on N. Finally, to complete the proof, we use the Hille-Yosida 
series expansion of the function M^ Up \t), t > 0. £ 

Remark 1. It can be shown that even linear SNNSS given by = Ao + d\k, /i*, = 
yU + d^k, k = 0, 1, . . . , s with d\ 7^ d^ have not the property stated in Proposition 2. 

In view of Proposition 2, we write w {l/p) (t) := w { ^ p \t), t > 0, < p < 1, N = s + 1, . . . , 
for the models (Ci) — (C4). The explicit expressions for the functions w^ Uv \t)^ t > is easy 
to obtain from (Ex)- (E 4 ). Now the Trotter - Kurtz approximation theorem ( ||14|| ) and (Ej)- 
(E 4 ) give immediately the following 

Corollary 2. For models (Ci) — (C 4 ) on a finite or infinite s-regular graph, 

Pr(^t v \x) = 1) = w M (t), t>0, xeV, 0<p<l. (4.103) 

In particular, in the presence of noise, the processes (Ci),(Cs), (C 4 ) have ergodic marginals, 
in the sense that for each of these processes 

\i mw M(t) (4.104) 

t— »oo 

exists and does not depend on < p < 1. 

The expressions (Di),(Ds) and (D 4 ) give correspondingly the following values for the limit 
in (gTUg) : 
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A + /i s ' 2 8h + a + b 

Remark 2. Formulae (D!)-(D 4 ) show the complicated influence of a constant additive noise 
on a transient behaviour of the process. In particular, note that under the absence of noise 
( h = 0) we have in ( |3.88| ), ( |3.91| ) a 2 = 0. It is also appropriate to mention that the processes 



(Ci)- (C4) are either attractive or anti-attractive. The latter means that in the definition of 
attractiveness ( [p4| , p. 132) the direction of inequalities for the flip rates is reversed. 

5 Ergodicity and Spectral gap 

Ergodicity. We will address the question of ergodicity of the processes (Ci)-(C 4 ) on a finite 
or infinite s-regular graph. The process (C 2 ) is, obviously, not ergodic. The processes (Ci), 
(C3) with a > and (C4) with a, b > are attractive. The key property of such processes is 
that ergodicity of their marginals implies the ergodicity of the process. A beautiful argument 
leading to this assertion is explained in |14[](see Corollary 2.8, p. 75, and Corollary 2.4, p. 136). 
So, by Corollary 2, we get 

Corollary 3. The following three processes are ergodic: (Ci) with hi + h 2 = A + [i s > 0, 
(C 3 ) with a > 0, h > and (C 4 ) with a, b > 0, h > 0. 

Remark 3 i.The ergodicity of the first among the three models in Corollary 3 was proven in 



| 12| . The ergodicity of the third one as well as the second one in the case s = 2, follows from 
the fact that these are attractive spin systems in one dimension with translation invariant 
and positive flip rates (see [|14]], Theorem 3.14, p. 152). To the best of our knowledge, the 
established ergodicity of the second model in the case s = 3 answers an open question. We 
will explain below that the e — M > condition in the case considered gives ergodicity for 
a < |/i only. 

ii. It is interesting to observe that, by Corollary 2, the models considered have ergodic 
marginals also in the case when they are not attractive and even not ergodic. For example, 
this is true for (C3), when h + a = 0, s = 2, in which case the process on Z\ is not ergodic, 
having two different absorbing states rji : r]\{x) = 0.5(1 + (— l) l+ ' x '), i — 1,2, x £ Z\. 
Spectral Gap. Recall (see for references 0], 0) that the spectral gap a > of an expo- 
nentially ergodic Feller-Markov process (ft t > with an invariant measure v on state space 
X = {rj} is defined by 

a = sup{[3 > : sup \Ef((pi v) ) - [ fdv\ < A f exp(-/%), / e C(X), t > 0}, (5.106) 
776* J 
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where Af > does not depend on t > 0. 

It is plain that formulae (Dj.), (D 3 ) and (D 4 ) for the mean coverage function provide an 
estimate from above for spectral gaps of the corresponding processes. Namely, if a is the rate 
of exponential convergence to the equilibrium of the function AfW(t), t > (=w("*)(t), t > 0), 
then a < a. From the other hand, the celebrated e — M > condition of ergodicity ( ||14|| , 
p. 31) provides the lower bound of the spectral gap of any ergodic IPS on a finite or infinite 
graph: a > e — M > 0. In the case of a SNNSS on a s-regular graph the quantities e and M 
are given by 

e = min {A fc + /ifc} M = s max {\X k - A fc -i|, \nk - fJ>k-i\} ■ (5.107) 

0<k<s 0<k<s 

Below are the values of e — M and a for the three ergodic models considered. 
(Ci): 

e-M = \ + n s = a (5.108) 

So, 

a = A + /is, (5.109) 
if Aq + yU s > 0. This fact was observed in |12| . Moreover, it was proven in [lOj that if G is a 



complete graph, then the noisy voter model is the only one SNNSS with the property (|5.109 ). 
(Cs): 

a = oti, (5.110) 

where ct\ is given by (|3.88|) , while 

{2h — sa, if a > . 

5.111 
2h + {s + l)a, ifa<0, 

where s — 2, 3. This says that in the case a > the e — M > condition is applicable for 

s — 2, if h > a, and for s — 3, if h > 1.5a. It is easy to verify that a± = e — M > iff a = 

which is the trivial case of the noisy voter model. 

(C 4 ): 

a = a u (5.112) 

where ot\ is given by ( |3.91| ), while 

e - M = 2h + min{0, a, 6} - 2{max \a\, \b\}, (5.113) 

where h = > 0. 
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Here again one can see that a,\ = e — M > is impossible unless a = b = 0. 
The preceding discussion leads us to the following 

Conjecture. Noisy Voter model is the only one NNSS for which a = e — M > 0. 

Concluding remark. It is demonstrated by our Theorem, that passing from the first order 
differential equation to the one of the second order does not enrich much the class of 
solvable ( in the sense of the mean coverage function ) SNNSS. A natural question arising 
in this connection is the characterization of SNNSS in the case of higher order differential 
equations. One can expect that the progress in this direction will lead to the discovery of 
wider classes of solvable SNNSS. The solution to this problem requires the analysis of the 
structure of generators gi, as defined in ( |3.6|) , of higher orders. At the moment the problem 
looks intractable. 
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